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Abstract 

Noncommutative features are introduced into a relativistic quantum field theory model of nu¬ 
clear matter, the quantum hadrodynamics-I nuclear model (QHD-I). It is shown that the nuclear 
matter equation of state (NMEoS) depends on the fundamental momentum scale, g, introduced by 
the phase-space noncommutativity (NC). Although it is found that NC geometry does not affect 
the nucleon fields up to 0{g 2 ), it affects the energy density, the pressure and other derivable quan¬ 
tities of the NMEoS, such as the nucleon effective mass. Under the conditions of saturation of the 
symmetric NM, the estimated value for the noncommutative parameter is rj = 0.014MeU/c. 
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I. INTRODUCTION 


Introducing noncommutative geometric features is believed to be an interesting way 
to generalize quantum mechanics |Mj. Noncommutative quantum mechanics (NCQM) 
arises as deformations of the Heisenberg-Weyl algebra. NCQM lives in a 2d-dimensional 
phase-space, where time is assumed to be a commutative parameter, and coordinate and 
momentum variables obey a NC algebra PE]: 

[x 1 , ad] = i6e l \ [p 1 ,^] = , [z^p*] = ih e ffS lj , (1) 

where e i3 is an antisymmetric matrix 1 , i and j range from 1 to 3, and 

n ‘fi = h ( J + il) - (2) 

The NC parameters, 9 and rj, are believed to be new fundamental constants of nature 
together with Planck’s constant. The fundamental scales of NC geometry are lg = \fd and 
l v = y/rj which must be obtained from experimental data for specific systems; for instance, 
the following bounds have been set 

lg < 2 x 10 _5 /m, l v < 8 x 10 _1 meU/c. (3) 

The NC extensions of QM show an impressive range of implications: on the quantum 
Hall effect [8], on the Landau level and the 2D harmonic oscillator problems in the phase- 
space pun!, and as a probe for quantum beating and missing information effects mi and as 
a source for quantum entanglement [12]. NCQM also admits violations of the Roberston- 
Schrodinger and Ozawa’s uncertainty relations [T3j. Furthermore, in the framework of 
quantum cosmology, phase-space noncommutativity has shown to give origin to the novel 
features for the black hole singularity mm, as well as to the equivalence principle [18] . 
One also expects some implications to compact objects [19]. 

In the present work, we examine the impact of NC in nuclear physics, in particular, 
on relativistic nuclear matter (NM) calculations. Nowadays, there is a growing interest in 
applications of the primary NM theories [2(11 [25] to study compact stellar objects. Of course, 
this has bearings on the mass-radius relationship and on studies of the crust thickness 
of neutron stars (NS), and on conditions for the collapse of NS and black holes. Other 

1 It is used that \zijktij = &k, £ijk^ij and e*, are Levi-Civita tensors and a normalized vector, respectively. 
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implications include supernovae explosions, the study of energetic heavy ions collisions 
[26], and on the properties of ordinary nuclei m 

In order to examine the NC effects into the NM calculations, we consider the QHD-I 
(or the era;-) model, a well known renormalizable relativistic quantum field model of the 
nucleon (p, n) system, interacting with the neutral scalar and vector mesons, cr, and, c o. 
Moreover, following Ref. [2%] . we assume that at high baryon densities, the scalar and 
the vector fields are replaced by their expectation values, which serve as a mean field in 
which the nucleons move. As will be shown, for a suitable value of rj under the empirical 
saturation conditions, the effective mass of a Dirac nucleon, M* (cf. Eqs. (20) and (37)), 
has a reasonable value. 

This paper is organized as follows. In section II, we briefly set the noncommutative 
tools to be used in our nuclear model. The field approach required to tackle the nuclear 
problem, the QHD-I model, is discussed in section III. In section IV we consider the NC 
QHD-I model. Finally, in section V, we present our results and discussion. 


II. SOME NC GEOMETRY TOOLS 


In order to generalize the QHD-I model subject to the algebra Eqs. (1), we must consider 
(see e.g., Ref. [29]) that NC fields satisfy the generalized Moyal * product 16): 


f(x)*g(x) = f(x) exp \ \ g(x) « f(x)g(x ) + ^0 ij d i f(x)d j g(x) + O(0 2 ) 


(4) 


This truncation will be sufficient for the purposes of our study and the dependence on the 
ij parameter will arise through the so called Seiberg-Witten map [30] . This is a noncanon- 
ical transformation, (x,p) i—> ( x',p '), that maps the NC algebra into the Heisenberg-Weyl 
algebra [6]: 


[x'ipx'j] = 0, [pfpjl = 0, Wi,Pj] = ittiij. (5) 

The NC variables, Eq. (1), can be mapped into the commutative ones, Eq. (5), through 
the Seiberg-Witten map: 

2 L) £ « p i' p ‘ = bp > + im) £ ‘V (6) 

where without loss of generality we choose a=b=l mm- 

In what follows we adapt the phase-space NC treatment of Dirac fields developed in Ref. 
[22] to the QHD-I model. As will seen this will allows us to directly assess the effect of the 
dependence on the NC parameter on the well known NM quantities. 
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III. A BRIEF REVIEW OF THE QHD-I MODEL 


We briefly review here the QHD-I ( au -) model (see also Ref. |28j). In the QHD-I model, 
the neutral scalar meson held, $, couples to the scalar density of baryons, TT, through 
the Yukawa interaction term with strength given by the coupling constant g S) while 

the neutral vector meson, couples to the conserved baryon current, 'I' 7 /i 'I', through 
g v ty with the coupling constant g v . In the mean held approach (MFA), baryons are 
assumed to move in a box of volume D within the mean held of the expectation values of 
the constant and condensed scalar and vector helds, $0 and Vo, respectively. The effective 
mass of a Dirac nucleon, M*, is given by M* = M — g s Q 0 , and the Lagrangian density for 
QHD-I model, in the MFA, is as follows [ 28] : 

Cqhd-i = V (i^ - Pg v V Q - M*) T + C 0 , (7) 

where Co is written via the constant scalar and vector meson mean helds [28] 

C 0 = ^m 2 ®,l - ^m 2 v Vj. (8) 

In Eq. (7), 7 m = ( 70 , - 7 ;), and 

7o = P, 7* = Pou- (9) 

In Eqs. (9), the Dirac matrices, P and cq, satisfy the following relations: 

[ai, dj] + = 2Sij, \ati, P] + = 0, a 2 = P 2 = 1, (10) 

so that 

0 a, 

CCj = 

<7 1 0 

where 1 and cq are the well known 2x2 unit and Pauli matrices, respectively. Using Eqs. 
(7) and ( 8 ), through the held equations, we obtain the following relations and the equation 
of motion for T: 


as 9 s 

$0 = ~^Ps, 

(12) 

m; 

T T 9V 

Vo — — 0 PBi 

(13) 

m; 

(if*^ - Pg v Vo - AT) T = 0, 

(14) 


1 0 
0 1 


( 11 ) 
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where p s = and ps = (4bT) are the scalar and baryon density, respectively, and 

the brackets <> denote the vacnnm expectation values. Since the fields are assumed to 
be constant, there is a static and uniform set of particles. Thus, considering free nucleons 
-0(k, A ){exp [fk.x — ie(k)t]} with momentum k, energy e(k) = e(|k|), polarization A, and a 
four-component Dirac spinor, -0(k, A), Eq. (14) leads to 

[ack + jdM*] ip (k, A) = e* ( k ) ip (k, A), (15) 

where e*{k) = [e[k) — g v \ o] is the effective energy of a nucleon. Regarding the superposition 
of solutions of Eq. (14), its general solution is given by: 

(*, t) = ~^= {^kA^(k, A)e* [k ' x - £+(fc)t] + Bl x V{ k, A) e -4k.x+ £ _(fc)t] j ^ ( 16 ) 

k,A 

where U(V) is the positive (negative)-energy spinor. A straightforward calculation, shows 
that: 


e ± (k) = [gM±E‘(k)], 


(17) 


where E*(k ) = \/k 2 + M* 2 . It should be pointed out that in Eq. (16), the summation 
over |k| is limited to kp, the Fermi momentum of the nucleons. On the other hand, 
the summation over A comprises summation over spin and iso-spin. Finally, through the 
familiar relationship between the energy-momentum tensor, energy density and pressure, 
one can write the NMEoS in the QHD-I framework 


£c = 


Pc = 


m v 


~PB 


1 

+ 2 


— (M 
9s 


1 2 


M* 


v 


r*k,F 


+ -^ I E*(k)k 2 dk, 


9v V 1 
- P B ) _ b 

m v / 2 


— (M-M*) 

9 s 


+ 


-(—) 

3 V2tt 2 / 




k 4 


E*(k) 


dk , 


(18) 


(19) 


where u, the degeneracy on the spin and iso-spin of nucleons, is 4. Of course, Eq. (19) 
follows from Eq. (18) through the relationship p = p\ \d[ejps) /dps]- The energy of an 
isolated system is made minimal at fixed volume, fl, baryon number, B, and temperature 
(here vanishing), by minimizing e with respect to M*. Now, M* is given by the self- 
consistency relation: 


M* = M - | — 
m x 


v 

2^2 


fkp l\/f* 

e^Lpdk. 


’0 


E*(k) 


( 20 ) 
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IV. THE NC QHD-I MODEL 


The NMEoS, Eqs. (18) and (19), as well as the nucleon effective mass, Eq. (20), are 
obtained from the Lagrangian density, Eq. (7), in the MFA approach. These relations are 
changed by the NC algebra, Eqs. (1). Following Ref. [29], the ordinary product of fields 
in the Lagrange density is replaced by the Moyal product, *, for the NC fields, T 7 . In the 
MFA, the Lagrangian density, in the NC geometry is as follows: 

C'qhd-i = ^ * (i^d, - Pg v V 0 - M*) * T 7 + C 0 . (21) 

As described in Section II, the NC fields and variables must be mapped into commutative 
ones. Thus, using the transformations Eq. (6), we get for the free noncommutative Dirac 
fields: 


T 7 (x 7 , t) ~ e 


ik'.x 


( 22 ) 


where it has been used that of k[ — p' t — —id\. Hence the NC Lagrangian density of QHD-I 
reads: 

C'nc = ~ \r) k H k x' 1 - (3g v V 0 - M*)T 7 + C 0 , (23) 

where mfd l = i^ l d[ — ^'rjkr/ k x n , and 6 = 0 was set in order to preserve gauge invariance 
(see Ref. [29]). Therefore, the equation of motion for the Dirac field reads 

- l -r] k a k x a - i3g v V 0 - M*\ = 0. (24) 

Eq. (15) is changed to 

cx.k 7 + -ck x x 7 .?7 + l3M* 

where s'* = e' — g v V q. For the spectral energy of the positive and negative states, e' ± , one 
finds: 


if'(k',\) = e'*(k')' l f'(k',\), (25) 


= 


g v V 0 ± E* 

with E'* = \J k 72 + M* 2 + k 7 x xbr/. For the energy density one finds: 


(26) 


£ - 2 {^ Pb ) + 2 


m 


— (M-M*) 

9s 


v 

+ SP 


E 

|k'|<fe F ' 


E,;dx 7 . 


(27) 
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In a similar way, one can get for the pressure: 


\ v / j |k'|<fc F v 


Expanding E'* up to 0(g 2 ), and converting the sum over k 7 into an integral, hence: 

-j / \ 2 -t - “2 

£ ' = o ( ~Pb ) +o —(M — M*) + —^ / E*(k')k' 2 dk' + 77 ^, (29) 


2 V m r 




2 i { v \ r kF k ,A 

+ 3V2t t 2 ) J 0 E*{k') 


where 


r 1 y 1 

£ 2(2t r) 3 fi 


k! x n dxdk ', 

E*(k') 


T = - — 
p 12 (2 tt) 3 n 


J,/3 7 / /k / 

[£*( A ;')] 3 ’ 


where E*{k') = vAc ' 2 + M * 2 and 0^ = (e^ x e x >).e v , where e*,/, etc. are the unit vectors 
in the direction of k 7 , etc. 

Choosing rj = 0, we recover Eqs. (18) and (19). The relevant quantities for the NC 
NMEoS are given by: 

£jvc = e C + V^et (33) 


PNC = PC + h r p, 


where 


Te = (Ef - 3M*E* F + 2M* 3 ) A, 


T„ = - ——^ - 2M* 2 ) E* f - - M* 2 — + -M* 3 1 A„, (36) 

p 12 (2vr) 3 [ V F ’ F 3 F Ep 3 J v ; 


12 (2tt) 


with EJj, = y / fc|T+ M * 2 and A v — f kl n x'dx 1 dEl k t/kl, where dily is the solid angle cl¬ 
ement of k!. The factor X v depends on the dimension of the box, here referred to as 
"noncommutative geometry length" ( NCGL ). 

The key quantity in the study of the NM properties is the nucleon effective mass. The 


effective mass in the NC case is given by: 


/ n \ 2 v f kp M* 

= iFM*' + " rM- ’ 
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where 


r m* = ( — 


i/ M* 



k'x' 


rdx'dk.'. 


^m s J (27t) 3 4Q J J [ E*(k')] 3 
As in Eqs. (33) and (34), the nucleon effective mass in the NC geometry is given by: 


Mfrc — Mq + r)T m* , 


where 




m « 


(2tt) 


-M 


*2 


E 


M 


. - .. -2 A,, 


M* 


E* 


(38) 


(39) 


(40) 


In the next section we present some numerical estimates for the NC effects. 


V. RESULTS AND DISCUSSION 


We consider as input the saturation point values of NM in the usual QHD-I framework 
. We show in Table I, the standard values for the pertinent parameters, where = 
, Cl = and M is the average of proton and neutron masses. 


Table I Relevant values for the NM in the QHD-I framework [28]. 


M(MeV) 

m v (MeV) 

m s (MeV) 


C 2 

938.93 

782.6 

550 

195.5 

267.1 


To investigate the effect of on the NM calculations, we consider the binding energy, 
£b = e/p — M, of the symmetric NM in which p p = p n = Pb/2- Considering the trivial 
constraint Pb = p P + Pn, we solve Eq. (29) and the self-consistent Eq. (37) with the 
assumption that \ v = 1. Fig. [l] shows the saturation curves of the symmetric NM for two 
different values of p with respect to the coupling constants from Table I. For p = 0 the 
saturation point takes place at kp = 1.42/m _1 with the value eb = —15.75 MeV [28]. In the 
case of p = 1, the saturation lies at kp = 1.41/m _1 with the value £b = —12.86MeU. For 
p 7 ^ 0 , although we cannot conclude that the saturation point is achieved, the qualitative 
behaviour of the curve is kept. 

Fig. | shows the behaviour of M*/M versus kp for p — 0, 1 for the parameters of Table 
I. It can be seen that M* is greater in the case of p 0 for 0.5 < kp < 3.5 fmr 1 . As the 
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effective mass controls the stiffness of the NMEoS, the NC geometry softens the EoS of 
QHD-I model for the relevant values of kp. 

To further investigate the effect of NC, we consider the following strategy. Applying 
the experimental saturation data of NM, £ 0 b = —15.86 MeV and po = 0.16/m -3 , and 
considering the trivial constraint on pp, we calculate the coupling constants g s and g v , and 
Ad* through Eqs. (29), (30) and (37), for different values of Pncgl = V^ry Table II shows 
the values of parameter 7]ncgl and the coupling constants which saturate the symmetric 
NM. The presented values of pncgl have been calculated in the natural units, h — c—1. 
As shown in Table II, the value of M*/Ad, computed by solving the self-consistent equation 
under the empirical saturation conditions and suitable value of Pncgl, lies in the acceptable 
interval, 0.7 < Ad*/M < 0.8 [32] . 


Table II The values of tjncgl, At*/M, g s < f > 0 , g v Vo, and the dimensionless coupling constants, which 
provide a fit for the experimental saturation data (eoi> = —15.86 MeV and po = 0.16/m -3 ). 


Vncgl 

Ad*/Ad 

g s MMe V) 

g v V 0 (MeV) 

c? 

C , 2 

4.87 x 10 -3 

0.78 

206.56 

127.18 

155.58 

91.20 

4.91 x 10 -3 

0.68 

300.46 

215.18 

223.09 

154.26 


The parameters in Table II show how the NMEoS changes as a function of the funda¬ 
mental momentum scale, g. With respect to the magnitude of coupling constants from 
Tables I and II, the vector repulsive and scalar attractive parts of the symmetric NM in the 
case of NC are smaller than the usual case. This fact can be verified through comparing 
values of parameters g s Q o and g v \f for the NC case, Table II, and those of QHD-I where 
~ 400MeH and g v V o ~ 330 MeV [SHJ. The NC geometry reduces the interaction mag¬ 
nitude of propagating nucleons in constant scalar and vector fields. As can be seen, this 
reduction is approximately similar for both attractive and repulsive parts of the interaction. 

We can now compute the value of 7] for the symmetric NM system. In order to do this, 
let us estimate the magnitude of the geometric parameter \ v = f (hd dklk'/kl, where 

— (g/v x e,j./).e, r If we assume e ^ = e z >, e v = e y i, and obtain dl v = 1; thus, A v will 
be the double-integral over the solid angle element of k' and the triple-integral over the 
magnitude of the position vector. Since the center of heavy nuclei is a typical example of 
NM, the root-mean-square radius ( R rms ) of heavy nuclei can measure the volume of our 
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box, fh A straightforward calculation yields X v ~ 70 fm for R rrns — b.bfm. Therefore, an 
estimated value of r/ is 6.96 x 10~ 5 in the natural units, or l v « 0.014MeI//c. It should be 
noted that we use the value of rj corresponding to M*/M = 0.78. 

We conclude, as a point of principle, that imposing the NC geometry in the QHD-I 
model, modifies the nuclear calculations and reduces the magnitude of nucleon-nucleon 
interaction selecting a suitable value of the NC geometry parameter, r /. On the other hand, 
it leads to a softer NMEoS than that of the usual case. Since the NC geometry softens the 
NMEoS, we can expect, for instance, that NC in the mass-radius neutron stars calculations 
might lead to a smaller maximum mass than the usual case. Of course, the present scheme 
can be extended to the other nuclear systems, such as neutron matter and neutron stars, 
by including the effect of other particles such as leptons and hyperons. This might have 
relevant implications for the understanding of nuclear matter under astrophysical conditions 
(cf. Ref. [19]). 
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